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ABSTRACT 


An  investigation  of  stresses,  arising  from  a  thermal  field,  in 
electronic  packages  comprised  of  a  ceramic  chip  attached  to  a 
composite  board  by  a  solder  joint  was  undertaken.  Two  types  of 
solder  joint  connections,  a  leaded  device  and  an  unleaded  device, 
were  analyzed.  The  stresses  result  from  the  mismatch  of  material 
properties  of  the  tri-assembly  system.  Procedures  are  shown  for  the 
analyses  of  the  two  types  of  solder  connections.  The  analyses 
showed  that  the  thermoelastic  stresses  for  the  stiff  unleaded 
connection  are  significantly  greater,  by  a  factor  of  20  for  a  the 
particular  application  undertaken  in  this  report,  than  the 
thermoelastic  stresses  in  the  flexible  leaded  connection. 
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Introduction 


Electronic  chip  packages  are  comprised  of  several  components  with 
different  material  properties  and  geometries.  One  of  the  most  basic 
electronic  packages  is  the  tri-assembly  system  consisting  of  a  chip 
attached  to  a  board  by  a  solder  connection.  During  the  fabrication 
of  this  simple  unit,  the  soldering  process  subjects  the  device  to 
a  thermal  field.  Since  the  three  components  of  the  device  have 
different  material  properties,  in  particular  their  Young’s  modulii, 
Poissons  ratios,  and  thermal  coefficients  of  expansion,  the  mismatch 
in  these  properties  result  in  thermoelastic  stresses.  In  addition  to 
the  thermal  field  arising  from  fabrication,  thermal  fields  and  hence 
thermoelastic  stresses,  result  from  the  operation  of  the  device  itself. 

Not  only  is  the  mismatch  in  material  properties  important  in  the 
thermoelastic  stresses  that  develop,  but  equally  important  is  the 
manner  in  which  the  solder  connection  is  made.  For  a  given  set 
of  chip,  board  and  solder  properties,  an  important  question  'what 
is  the  best  way  to  make  the  solder  connection?'.  With  regard  to 
geometric  configuration,  there  are  two  main  types  of  connections, 
the  leaded  connection  and  the  unleaded  connection.  In  the  leaded 
connection,  the  chip  is  attached  to  the  board  by  a  thin  ’lead’  as 
shown  in  Figure  1.  In  this  configuration,  the  chip  is  somewhat 
removed  from  the  board  its  attached  to.  In  contrast,  in  the 
unleaded  connection,  the  chip  is  attached  along  its  edges  to  the 
board  substrate  by  a  thin  layer  of  solder.  This  unleaded  joint 
arrangement  is  significantly  ’stiffer’  than  the  leaded  connection,  and 
all  other  things  being  equal,  one  would  expect  that  the  stiffer 
unleaded  joint  would  result  in  larger  thermoelastic  stresses.  This 
investigation  seeks  to  determine  the  relative  magnitudes  of  the 
stresses  from  these  two  kinds  of  solder  joints.  It  should  be  noted 
however,  that  the  proximity  of  the  chip  to  the  board,  in  the 
unleaded  joint,  might  result  in  a  more  uniform  thermal  field  thereby 
reducing  the  thermal  gradients  and  thus  the  thermal  stresse*'^  The 
question  of  the  thermal  field  that  results  from  each  of  th**  solder 
joints  is  not  addressed  in  this  investigation  which  only  jonsiders 
the  effect  of  the  same  uniform  thermal  environment  o  i  each  of 
these  (leaded  and  unleaded)  solder  joints. 
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Analysis  of  leaded  device 

Figure  1  shows  the  geo'-netry  and  dimensions  of  the  leaded 
device,  considered  in  this  study  which  consists  of  a  ceramic  chip, 
lead  wires  and  a  circuit  board.  The  cross-sectional  dimension  of  the 
lead  wire  is  given  m  Figure  2.  » 

Three  general  assumptions  are  made.  First,  the  whole  device  is 
under  a  thermal  cycling  load  of  +65"c  ~  -20  c  and  the  lead  wire  has 
an  initial  strain  from  soldering  process.  The  device  is  soldered  to 
the  circuit  board  at  a  temperature  of  -i-185=c  and  this  tempature  is 
assumed  to  propagate  only  to  the  lead  wire  during  the  soldering. 

Another  assumption  is  on  the  deformation  of  the  device  that  the 
ceramic  chip  and  the  circuit  board  are  so  much  stiffer  than  the  leads 
so  all  the  deflections  due  to  the  thermal  loading  occur  only  to  the 
leads.  Finally,  the  device  is  assumed  to  be  geometrically  symmetric 
about  the  center  so  the  following  calculation  is  performed  only  for 
the  half  of  the  device. 

Material  properties  of  the  chip,  lead  and  the  glass-epoxy 
substrate  of  the  circuit  board  are  as  follows; 

Thermal  expansion  factor  of  the  chip,  =  6.5  lO'®  ppm/  c 

Thermal  expansion  factor  of  the  lead,  =  4.47  lO'^  ppm/x 

Thermal  expansion  factor  of  the  substrate,  =  16.  10-®  ppm/^c 

Young's  modulus  of  the  lead  wire,  E  =  21.  10^  psi 

Moments  of  inertia  of  the  lead  wire,  I  =  bh^/12  =  0.017  (0.005)^/12 

=  1.77  10-10  in'^ 


Calculation  of  thermal  mismatch 


At  -1-20  c; 


Initial  deformation  of  a  lead  wire  after  the  welding 
(from  -(-1 85  c  to  +20  c) 

=  a  AT(L)  =  4.47  10-0(-165  c)(0.045)  =  -33.2  lO'O  in 

s 
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At  +65‘c: 

Deformations  due  to  temperature  changes  from  +20'c  to 
+65‘c 

Achip  =  aj.AT(L^)  =  6.5  10-6(45=c)(0.26/2)  =  38.0  10-6  in 
A,ead  =aLAT(LL)  =  4.47  1 0-6(45“c)(0.045)  =  9.1  10-6  in 
Asub  10-6(45”c)(0.35/2)  =  126.0  10-6  in 

Mismatch  between  the  top(chip  and  lead)  and  the 
bottom  portion(substrate)  is 

^65^  ~  ^top  '  ^bot 

^■^chip  ^lead  ^INIT^  *  ("^sub^ 

=  (38.0  10-6  +  9.1  10-6  -  33  2  io-6)  -  (126.0  10-6) 

=  -112.1  10-6  in 


At  -20'c; 

Deformations  due  to  temperature  changes  from  -t-20"c  to 
-20=0 


A^hjp  =  acAT(L^)  =  6.5  10-6(-40=c)(0.26/2)  =  -33.8  10-6  in 

A|g3d  =  =  4.47  10-6(-40=c)(0.045)  =  -8.05  10-6  in 

Ag^P  =a3AT(L3)  =  16.  10-6(-40=c)(0.35/2)  =  -112.  10-6  in 

Mismatch  between  the  top(chip  and  lead)  and  the 
bottom  portion(substrate)  is 

*^-2C’  “  -^top  ■  -^bot 

^•^chip  -^lead  ■^INiT^  *  ^‘^sub^ 

=  (-33.8  10-6  -8.1  10-6  -  33.2  10-6)  -  (-1  12.  10-6) 
=  36.9  10-6  in 

Since  Aggc  is  much  greater  than  A_20'  '^s  magnitude,  it 

is  considered  that  Agg.  becomes  critical  for  the  stress 
calculation. 


AnalPQ.v.,  Qf.jLemDerature  change  effect  and  support  settlement  effect 
It  can  be  seen  tnat  tne  force  balance  from  the  temperature 
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changes  are  the  same  as  that  from  the  support  settlement  shown  in 
Figure  3.  So  we  now  treat  the  displacement  Agg,  from  the  above 
calculation  as  the  settlement  5. 


^  -  ^65= 

To  find  the  force  P  cause  by  a  displacement  5  which  is  known 
from  the  previous  calculation,  we'll  break  the  frame  ABC  into  beam 
AB  and  beam  BC  as  in  Figure  4  and  5. 

Denoting  and  as  the  length  of  the  beam  AB  and  BC, 
respectively,  the  tip  displacements  and  the  tip  rotations  at  the  beam 
AB  due  to  an  applied  force  Pg  and  a  moment  Mb  are  given  by  the 

formulas 


PL  ^ 

3EI 

V 

2EI 

(1) 

PL^^ 

ML,, 

=  2EI 

El 

(2) 

and  for  beam  BC  the  tip  rotations  can  be  expressed  by 

El 


(3) 


Compatibility  requires  that  0g  at  the  above  equations  (2)  and 

(3)  be  the  same  so  by  equating  these  two  equations  an  expression  for 
M  can  be  found. 


+  L,) 


(4) 


Substituting  (4)  into  (1) 


5  = 


L,  +  41^ 
12(L,  +  L,) 


or  P  can  be  written  in  terms  of  6. 
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*-v 


112.1  10-6  (21.  106)  (1.77  10-6)  12(0.028  +  0.045) 
(0.028)^  0.028  +  4(0.045) 


=  0.0811  lb 


Shear  stress  calculation 


The  shear  stresses  at  the  solder  joint  are  obtained  by  dividing  the  force 
by  footage  area  of  the  lead  wire  which  is  soldered  to  the  circuit 
board, 

P  _P  0-0811 

^  ■  A  "  bL  ~  0.017(0.03)  ■  P®' 
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Analysis  of  unleaded  device 

Figure  6  shows  the  geometry  and  dimensions  of  the  unleaded 
device  which  consists  of  a  ceramic  chip  fixed  to  a  circuit  board 
by  a  thin  strip  of  solder.  In  the  present  application,  the  solder 
adhesive  layer  covers  only  a  small  area  along  the  edge  of  the 
chip. 

Two  models  of  the  problem  follows.  The  first  model,  developed 
by  Suhir  (1),  provides  an  analytical  determination  of  thermally- 
induced  stresses  in  a  tri-material  assembly.  The  second  model  is 
a  finite  element  mode!  which  utilizes  a  recently  developed  element 
which  maintains  inter-layer  continuity  along  the  chip-solder  and 
board-solder  interfaces. 

Suhir's  Model 


A  free  body  diagram  of  the  tri-material  assembly,  showing  the 
notation  used  hereafter,  is  shown  in  Figure  7.  The  Suhir  model 
is  applicable  to  assemblies  with  continuous  attachment  (a  =  I), 
as  well  as  non-continuous  attachment  (a  <  I).  The  development 
is  carried  through  by  invoking  equilibrium,  displacement 
compatibility,  and  thermoelastic  constitutive  relations.  The 
displacement  compatibility  relations  impose  the  condition  of  inter¬ 
layer  displacement  continuity.  For  a  detailed  development  of  the 
analytical  relations,  the  reader  is  referred  to  the  original  work. 
Here  a  brief  summary  of  the  development  is  presented,  followed 
by  a  sequential  list  of  the  calculations  that  lead  to  the  maximum 
interlayer  shearing  stress. 

Denoting  the  axial  displacements  at  the  top  and  bottom  of  the 
i'^  layer  as  u*,(x)  =  u  ,(x)  respectively,  the  displacement 
compatibility  conditions  along  the  interfaces  of  layers  1-2  and 
layers  2-3  are: 

ul  (x)  =  U2  (x)  (5) 


ui  ix)  =  U3*(x)  (*) 

Axial  displacements  u,(x),  u*2(x),  U2(x),  and  u''3(x)  are  obtained 
by  adding  the  individual  contributions  to  these  displacements 
arising  from  (a)  thermal  expansion,  (b)  axial  forces,  (c)  shear 
forces,  and  (d)  bending  moments.  The  moments  are  related  to 
the  axial  forces  by  the  moment  equilibrium  equation  which  gives 
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(7) 


2 


T^(x} 


h^+h 

~ 


^  Tj(x)  =  *  M,  * 


where  T,  and  Tj  are  the  axial  forces  on  layers  1  and  2.  The  T| 
axial  forces  are  related  through  equilibrium  to  the  shear  stresses 
between  layers,  that  is 


T,{x)  =  dC 


(8) 


where  t,(x)  is  the  shear  stress  between  the  i"’  and  (i  +  1)*‘ 

layers.  The  M,  moments  are  related  to  the  radii  of  curvature,  p, 
through 

M,  (X)  =  -  — ^  i  =  1,2,3  (9) 

p  (x) 


and  the  C,  flexural  rigidity  coefficients  for  each  layer  are  given 
by 


where  E,,  h,,  and  v,  are  the  Young’s  modulus,  thickness,  and 
Poisson's  ratio  of  the  i'^  layer  respectively.  Equation  (10)  assumes 
there  is  no  variation  of  lateral  displacement  with  respect  to  y, 
and  therefore  the  radii  of  curvature  of  layers  are  all  equal. 


Using  equilibrium  equation  (7)  and  moment-curvature  equation  (9) 
yields 

-  T  (x)  (11) 

where 


D  =  D, 


D, 


(12) 


is  the  flexural  rigidity  of  the  assembly,  and  where 


.1*1 


2D 


(13) 


IS  an  additional  compliance  due  to  bowing.  The  individual 
contributions  to  Eqn's  (5)  and  (6)  are  the 
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•  the  thermal  contribution  is; 


(AT) 


(14) 


•  the  axial  force  contribution  is: 

-xJ"T,(^)di 

0 


(15) 


where  /,  is  the  axial  compliance  of  the  i  layer,  given  by 


X  =  _ 'ii 

^  E.h, 


(16) 


•  the  shear  stress  contribution  is:  -A:,r,(x),  where  a:,  is  the 
interfacial  compliance  given  by 


K. 


h. 

3^ 


(17) 


•  the  rotational  contribution  is: 

2  Jo  p(0 


(18) 


Invoking  the  condition  that  the  shear  stresses  be  anti-symmetric 
with  respect  to  the  origin  of  coordinates  leads  to  the  expression 

Tj  (x)  =  C,  sinh  icx  i  =  l,2 


where  k  is  obtained  by  solving  a  characteristic  equation  defined 
in  Eq.(20)  below.  Substitution  of  the  above  expressions  into 
compatibility  conditions  (5)  and  (6)  yields  the  quartic  equation 

K..K,,k*  -  +  ft  =  0  (20) 

where  the  k^^'s  and  A.^'s  are  interfacial  and  axial  compliances, 
respectively,  of  the  tri-material  assembly,  given  by 

K.^  =  K.  +  <2  (21) 
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•^23  "-s 

(22) 

>-°23  =  ^2  ^3  '*■  ^23 

(23) 

A,i2  ~  ^  ^2  ^  ^12 

(24) 

6  =  X3A.2  +  ^2^3  ^3^1  *  ^3^12  ^2  ^^12  *  V^^23  ^  ^ 

(25) 

Letting 

0)  =  k^,  Eq.(20)  becomes 

-  (K.-A.23  +  K23X°2)t^  +5  =  0 

(26) 

whose 

solution  is 

-b±\Jh'^-4ac 

2  ~ 

^:a 

(27) 

where 

a  =  K.jX,,  ^  =  -(K32>.°3  '  K23XI2) 

(28) 

The  k, 

(i  =  1,2)  roots  are  then  obtained  as 

^1.2  =  ^'^3.4  =  ±V''^ 

(29) 

The  expressions  presented  above,  valid  for  any  tri-material 
assembly,  reduce  to  the  simpler  form 

T  (x)  =  k-  ^^~-J'-  sinhkx  (30) 

Acoshi 


where 


K  =  +  2X3  +  K3 


(31) 


when  the  Youngs  modulus  of  the  mid  layer,  or  its  thickness  is 
very  small,  as  in  the  case  of  a  thin  layer  of  solder  used  as  an 
adhesive.  In  accordance  with  Eq.(19),  the  maximum  shear  stress 
occurs  at  X  =  f,  and  can  be  calculated  from 
where 
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(32) 


X 


max 


1 


cosh  k{  1  -  a) 
cosh  kl 


(33) 


Application  of  Suhir  Model  to  Unleaded  Device 

Here,  the  closed  form  equations  from  Suhir’s  model  are  used  to 
determine  the  shear  stresses  in  an  uleade^  device.  The 
application  is  to  unleaded  device  with  dimensions  and  material 
properties  as  follows: 


Material  #1;  Ceramic  chip; 

fi,  =  3.6x10®  in/in/F,  v  =  0.3,  E,  =  40x10®  psi, 

G,  =  15.4x10®  psi,  h,  =  .056  in. 

Material  #2;  Solder: 

=  25.X10®  in/in/F,  v  =  0.3,  =  1.86x10®  psi, 

Gj  =  .66x10®  psi,  hj  =  .002  in. 

Material  #3;  (board): 

«3  =  16.X10®  in/in/F,  v  =  0.4,  E3  =  2.0x10®  psi 

G3  =  .77x10®  psi,  hj  =  .056  in. 

The  above  parameters  lead  to  the  following  compliances: 

Chip:  D,  =  643  lb. in..  A,  =  .3125x10®  in. /lb., 

x,  =  1.212x10®  in. /lb. 

Solder:  Dj  =  .0231  lb. in.,  Aj  =  64.52x10®  in. /lb., 

Xj  =  2.525x10®  in. /lb. 

Board:  D3  =  32.2  lb. in.,  Aj  =  6.25x10®  in. /lb., 

X3  =  24.24x10®  in. /lb. 

The  other  parameter  values  necessary  for  a  solution  are: 


A, 2  =  1.378x10®  in. /lb. 


Ajj  =  1.378x10®  in. /lb. 
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AT ,  2 

=  3.737x10® 

in. /lb. 

^2.3  — 

26.77x10®  in. /lb. 

A, 2' 

-  64.83x10® 

in. /lb. 

'^■23  ~ 

70.77x10®  in. /lb. 

=  66.21x10® 

in. /lb. 

^^23°  = 

72.15x10®  in. /lb 

Substitution 

of  these  parameters  into 

Eq.(25) 

yields 

d  = 

.7897x10® 

(in. /lb.)® 

which  upon 

substitution 

into  Eq.(26) 

yields  k. 

=  19.85  in.  V 

Finally,  the  maximum  shear  stress  at  the  solder  interface  is 

obtained  from  Eq.(32)  as 

=  4360  psi 

The  simpler  approximate  expression  for  k  from  Eq.(31)  yields  k 
=  19.61  in.'\  which,  in  turn,  gives  =  4310  psi. 

Finite  Element  Model 

An  alternate  model  of  the  tri-assembly  problem  using  a  finite 
element  method  was  undertaken.  Continuity  along  the  interfaces 
of  adjacent  layers  is  maintained  by  selection  of  the  six  degree 
of  freedom  element  shown  in  Figure  8.  The  six  degrees  of 

freedom  are  as  follows: 

(),  is  the  axial  displacement  at  the  lower  left  corner 

is  the  axial  displacement  at  the  upper  right  corner 

(^3  is  the  lateral  displacement  at  the  left  end 

is  the  axial  displacement  at  the  lower  right  corner 

is  the  axial  displacement  at  the  upper  right  corner 

f’)s  is  the  lateral  displacement  at  the  right  end 

The  axial  displacement  field  u(x,y)  is  taken  as  linear  interpolation 
in  both  the  axial,  x,  and  transverse,  y,  directions.  That  is, 

u(x,  y)  =  (y)  (N.^  (x)  5^  (x)  82)  *  iy)  (x)  5^  (x)  b^) 
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and  the  lateral  (i.e. .transverse)  displacement  field  is 
function  of  x,  that  is, 

v{x)  -  Kix)  63  N^ix) 


where  the  linear  interpolation  functions  are: 

(y)  =  \ 


H-Jy)  -  ^ 


A’,  (x)  =  ^  (  i  -  x) 


h\(x)  =  ^ 


The  strain-displacement  relations  are: 


=  £l^  (6,-63)  (6,-6,: 


<9x 


(3u  N'.(x)  N^{x)  6,-6. 

Y^.  =  3-  ■^3-  =  —  (6^-6,)  +  ■ - 

-  dy  dx  n  ^  h  -  ^  1 


On  the  element  level,  we  have  the  bending  effects 

^  B  ~ 

where 

6s  =  <63  6,  6,  65> 

and 

-  j  ^  j  ^BEB'^dydx 


linear 

(35) 

(36) 

(37) 

(38) 

(39) 

(40) 

(41) 

(42) 

(43) 

(44) 
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and 
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The  bending  and  shear  element  stiffness  matrices  are  combined 
to  form  the  element  stiffness  matrix,  which  in  turn  are  again 
combined  to  form  the  system  stiffness  matrix.  In  accordance  with 
standard  FEM  techniques,  the  system  matrix  equations  are  solved 
for  the  system  displacement  degrees  of  freedom.  The  element 
strains  are  then  calculated  from  Eq.(40)  and  (41)  above.  Finally, 
the  average  bending  stresses  at  the  bottom  and  top  of  an 
element  are  obtained  from  the  stress-strain  relations 


6,-6, 


6^-6. 


Similarly,  the  average  element  shear  stresses  are  calculated  from 
the  strress-strain  relations 


6,-6.  6,-6, 

^  *  ■  - 

h  1 


)  =  J  ^ 

1/7  “T 
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The  finite  element  mesh,  shown  in  Figure  9,  consisting  of  104 
elements  and  149  degrees  of  freedom  produced  an  interface 
maximum  shear  stress  of  3650  psi.  This  FEM  calculated  stress  is 
at  the  midpoint  of  the  right  most  element.  This  stress,  obtained 
from  an  FEM  analysis  is  at  the  midpoint  of  the  interface  at  the 
right  most  element.  Quadratic  extrapolation  to  approximate  the 
hyperbolic  sine  expression  for  the  r  shear  stress  brings  the  FEM 
predicted  value  to  about  3800  psi.  Thus  the  predicted  values  of 
the  maximum  shear  stress  from  the  Suhir  model  and  the  FEM 
model  differ  by  less  than  15  percent. 
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Conclusions 


This  investigation  has  shown  that  under  the  same  thermal 
environment,  the  stiff  unleaded  solder  joint  has  significantly  higher 
thermoelastic  stresses  than  the  flexible  leaded  solder  joint.  Since 
the  unleaded  joint  has  a  stronger  suppression  of  displacement 
than  the  leaded  joint,  and  thermal  stresses  result  from 
displacement  constraint,  the  result  is  not  an  unexpected  one.  In 
the  particular  application  presented  in  the  previous  sections,  the 
stresses  in  the  unleaded  connection  were  twenty  times  the 
stresses  in  the  leaded  connection.  This  large  degree  of  difference 
between  the  stresses  in  the  two  cases  is  surprising. 

As  previously  noted,  the  comparison  is  for  the  case  where  both 
connections  are  subjected  to  the  same  thermal  environment.  This 
assumption  does  not  take  into  account  the  fact  that  the  thermal 
field  is  very  likely  to  be  different  for  the  two  types  of  joints. 
One  might  expect  that  the  leaded  joint  with  the  'separated' 
components  may  have  larger  thermal  gradients  than  the  unleaded 
joint  where  the  components  are  in  close  proximity  to  each  other 
and  therefore  may  have  a  more  uniform  thermal  field.  In  fact,  the 
solution  to  the  stress  problem  really  requires  its  coupling  to  the 
thermal  (energy  balance)  problem.  This  later  coupled  (and 
transient)  problem  is  presently  under  investigation  as  a  NPS 
thesis  topic  under  the  direction  of  D.  Salinas  and  Y.W.  Kwon  (2). 


16 


substrate 


.03  in 


Figure  1,  Geometry  of  leaded  device 


Figure  2.  Cross  section  of  lead  wire 
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Figure  3.  Support  movement  <*► 
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Figure  4,  Deformation  of  beam  AB 


Figure  5.  Deformation  of  beam  BC 
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Figure  7.  Force  and  moment  equilibrium  of  components 


2 


References 


1.  E.  Suhir,  ’Chapter  5,  Thermal  Stress  Failures  in  Microelectronic 
Components  -  Review  and  Extension',  pp. 337-412,  in  Advances  in 
Thermal  Modeling  of  Electronic  Components  and  Systems,  Volume 

1,  edited  by  A.  Bar-Cohen  and  A.D.  Krauss,  Hemisphere 

Publishing,  Corp.,  New  York,  1988. 

2,  D.  Salinas  and  Y.W.  Kwon,  ’Thermoelastic  Stresses  in  Multi¬ 
layered  Media  in  a  Non-uniform  Temperature  Field’,  submitted  for 
presentation  to  the  Pan  American  Congress  of  Applied  Mechanics 
(PACAM  III)  to  be  held  in  Sao  Paulo,  Brazil,  Jan  4-8,  1992. 


DISTRIBUTION  LIST 


Office  of  Naval  Research 
Director,  Technology  (Code  02) 

800  North  Quincy  Street 
Arlington,  VA  22217-5000 

Office  of  the  L>irector  of  Defense 
Research  and  Engineering 
Information  Office  Library  Branch 
The  Pentagon,  Rm.  5E  1006 

Defense  Technical  Information  Center 
Cameron  Station 
Alexanderia,  VA  22314 

Naval  Postgraduate  School 
Technical  Library  (Code  52) 

Monterey,  CA  93943-5100 

Research  Administration  (Code  81) 

Naval  Postgraduate  School 
Monterey,  CA  93943-5100 

Professor  A.J.  Healey 
Dept,  of  Mech.  Engn'g. 

Naval  Postgraduate  School 
Monterey,  CA  93943-5100 

Prof.  Young  W.  Kwon 
Dept,  of  Mech.  Engng. 

Naval  Postgraduate  School 
Monterey,  CA  93943-5100 

Commanding  Officer 

Naval  Weapons  Support  Center  (Code  0211) 
Crane,  Indiana  47522 

Prof.  Young  S.  Shin 
Dept,  of  Mech.  Engn'g. 

Naval  Postgraduate  School 
Monterey,  CA  93943-5100 

Prof.  Philip  Y.  Shin 
Dept,  of  Mech.  Engn'g. 

Naval  Postgraduate  School 
Monterey,  CA  93943-5100 

Prof.  David  Salinas 
Dept,  of  Mech.  Engn'g. 

Naval  Postgraduate  School 
Monterey,  CA  93943-5100 


1  copy  , 

I 

1  copy 

2  copies 

2  copies 

1  copy 

1  copy 

1  copy 

3  copies 

1  copy 

5  copies 

5  copies 

I 


22 


